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ABSTRACT 

We analyze the Minkowski functionals with a large A^-body simulation of a standard ACDM model, focus- 
ing on transition scales between linear and non-linear gravitational evolution. We numerically calculate the 
Minkowski functionals with sufficient accuracies to investigate the transition scales, 10-50/i~'Mpc. The re- 
sults are compared with analytic formulae of linear and second-order perturbation theories. We first show that 
the skewness parameters of the density fluctuations, which are important in second-order analytic formulae of 
the Minkowski functionals, are in good agreement with the perturbation theory. Considering relative differ- 
ences between the Minkowski functionals of the analytic formulae and that of the simulation data, we evaluate 
accuracy levels of the predictions of the perturbation theory. When the straightforward threshold i' by density 
value is used in Minkowski functionals, the accuracy of the second-order perturbation theory is within 10% 
for smoothing length R > 15/i~'Mpc, and within a several % for R > 20/i~'Mpc. The accuracies of the linear 
theory are 2-5 times worse than that. When the rescaled threshold by volume fraction, ff is used, accuracies of 
both linear and second-order theories are within a few % on all scales of lOr'Mpc <R< SOr'Mpc. 
Subject headings: cosmology: theory — large-scale structure of universe — methods: statistical 



1. INTRODUCTION 

The morphological patterns of the large-scale structure in the universe are the consequence of the dynamical evolution of the 
universe, and thus provide important information on cosmology. The pattern of the clustering of galaxies is described in many 
ways. For example, Soneira & Peebles ( 1978) indicated the picture of the hierarchical clustering, and Joeveer & Einasto ( 1986') 
indicated that of the cell structure. Gott, Melott, & Dickinson (1986) (hereafter GMD) proposed a new picture, which is depicted 
by sponge-like topology of the universe. To capture the topology of the universe, it is useful to look at the three dimensional 
views of galaxy distributions (Einasto & Miller 1983). The clustering pattern is so complicated that it is difficult to fully describe 
the pattern by either the hierarchical picture or the cell picture alone. 

GMD showed that the large-scale structure of the universe has sponge-like topology by analyzing the CfAl catalog of galaxies. 
They examined various cosmological models of the galaxy distribution and quantified the topology by the genus statistic. The 
CfAl catalog which GMD used is the redshift survey of about 1,800 galaxies, which is relatively small catalog with respect to 
the today's standard. After that survey, numerous redshift surveys such as the Las Campanas Redshift Survey (LCRS), the IRAS 
Point Source Catalog Redshift Survey (PSCz), etc. have been carried out. Recently, the 2-degree Field Galaxy Redshift Survey 
(2dFGRS) was completed. The Sloan Digital Sky Survey (SDSS) is now in progress. In the SDSS, redshifts of about 8 x lO'' 
galaxies are being observed. These surveys are large enough to quantitatively discuss the topology of the large-scale structure on 
large scales. Topological analysis is one of the major methods to analyze these redshift surveys (see, e.g., Vogeley et al. 1994; 
Colley 1997; Canavezes et al. 1998; Colley et al. 2000). 

The Minkowski Functionals (MFs) are also used as descriptors of the morphology of the large-scale structure. While they are 
originally mathematical quantities, Mecke, Buchert, & Wagner (1994) applied the MFs to the analysis of the galaxy clustering. 
Schmalzing & Buchert (1997) provided a computational algorithm to calculate the MFs. The genus statistic is one of the MFs in 
a certain condition. 

The MFs possess complementary information to the two-point correlation function or the power spectrum, which are the most 
fundamental tools to quantify the clustering pattern. One of the important applications of the genus statistic and the MFs is a 
Gaussianity test of the primordial density field, which can not be performed by the two-point correlation function or the power 
spectrum. Any Gaussian field has a universal form of the genus statistic and the MFs as functions of the density threshold. As a 
result, any deviation from that form indicates the non-Gaussianity of a density field. 

Hikaae et al. (2003a) analyzed the early SDSS sample using the MFs. They found the MFs of the SDSS are consistent with 
that of the ACDM model. However, their analysis is limited by the cosmic variance because of the smallness of the data, and 
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the discrimination of models with various parameters is still difficult. Complete SDSS data would distinguish different models in 
more detail. 

When we compare observations with theories, the A^-body simulation plays an important role. A cosmological A^-body sim- 
ulation was first carried out by Miyoshi & Kihara (1975). Since then, the number of particles in A^-body simulations has been 
growing every year Nowadays simulations with the number of particles of over 10^ and the physical size of 1 /i"'Gpc^ have come 
to be available. Springel et al. (1998) compared the genus calculated from an A^-body simulation, which contains 256^ particles 
in a (240/i~'Mpc) box, with the data of the IRAS 1.2-Jy redshift survey. Hikage et al. (2002) used an A^-body simulation of the 
ACDM model carried out by Jing & Suto (1998), which contained 256^ particles in a (300/;~'Mpc)^ box. They found that the 
genus statistic of their A^-body simulations is consistent with the one calculated from the data of the SDSS Early Data Release. 
Recentlv. lHikage. Taruva. & Su to (2003b) used the Hubble volume simulation with a box size of (3000/i-'Mpc)3 and a number 
of particles A^ = 10^, and investigate the genus statistics with respect to the biasing effects. 

A drawback of using the A^-body simulation in comparing theories and observations is that the A^-body simulation is compu- 
tationally costly to obtain theoretical predictions of given cosmological models. Fortunately, there are analytic approximations 
for MFs. The MPs of a random Gaussian field is well-known (Tomita 1986). However, the density field in the universe is not 
exactly random Gaussian even if the primordial density field is Gaussian, because of the nonlinear evolution effects. Matsubara 
(2003) derived the analytic approximations of the MFs of weakly non-Gaussian fields, and obtained analytic formulae of the MFs 
with effects of the weakly nonlinear evolution by applying the second-order perturbation theory. Since in the last formulae is 
assumed an approximation that the non-Gaussianity is weak, there is a regime that the formulae can be applied. Identifying such 
regime is indispensable to use the analytic formulae in analyses of the observations, instead of performing computationally costly 
simulations model by model. 

In this paper, we analyze a large A^-body simulation of the large-scale structure using the MFs and investigate the transition 
scales between linear and non-linear evolution, and compare them with the theoretical formulae. The accuracy levels of the above 
analytic formulae are identified. We use an A^-body simulation of 512-^ particles in a box of (1024/!"'Mpc)-'. 

This paper is organized as follows. In §2, we briefly review the concept of MFs and their theoretical formulae. In §3, we 
explain the computational methods of the MFs. In §4, we analyze the MFs of the simulation data and compare them with the 
analytic formulae. In this section the comparison of the skewness parameters, which play an important role in weakly nonlinear 
formulae, is also presented. In §5, we summarize the results and discuss future outlook. 

2. THE MINKOWSKI FUNCTIONALS 
2.1. Definitions 

Mecke, Buchert, & Wagner (1994) and Schmalzing, & Buchert (1997) proposed the MFs as the geometrical descriptors of the 
distribution of galaxies. These quantities are originally derived by Minkowski (1903). To define the MFs of galaxy distributions, 
we first consider isodensity contours of the density contrast 5 = (p- p)/ p, where p is the density field and p is the mean density. 
Then we identify the set of regions M where the density contrast 5 exceeds some threshold. We use the threshold v which 
is defined by = (^th/co for the moment, where 5(a is the threshold density contrast and do = \/ {5-) is the rms of the density 
contrast. There are t/ + 1 MFs in a J-dimensional space. The case li = 3 is of our primary interest. In this case, the MFs correspond 
to the following quantities: (1) the fractional volume enclosed by contours 

Vo=l- j dV, (1) 

where V is the total volume, (2) the surface area per volume 

^1 = ^ / "^'^^ 

where the region of the integral dM is the surfaces of the contours, (3) the integrated mean curvature 
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where 1 /Ri and l/Ro are the principal curvatures on the surface, and (4) the integral Gaussian curvature 



1 f 1 



^3 = T^/ ^T^d^A. (4) 
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The last quantity is proportional to the Euler characteristic of the contour surfaces. The Euler characteristic is also propor- 
tional to the genus statistic if the boundary of the sample is neglected. Further geometrical meanings of MFs are found in 
iMecke. Buchert & WagneJ ( 11994 . 

2.2. Prediction for First- and Second-order Perturbations 
2.2.1. Gaussian random field 

The analytic formulae of the MFs as functions of the threshold v for random Gaussian field is given by (Tomita 1986; Schmalz- 
ing & Buchert 1997); 

1 / 1 \ 

V.(^)= . (p.) H,.,iv)e--I\ (5) 



where k = Q, 1,2,3, and the factor uik is the volume of the unit ball in /:-dimensions, i.e., luq = 1, wi = 2, uj2 = tt, and = 4tt/3. 
The quantities ctq and aj are defined by 
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and the functions H„{n = 0,2,3, • • •) are the Hermite polynomials with a convention, 
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Following Matsubara (2003), we use a notation 
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when k = 0. 
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2.2.2. Non-Gaussian field 

In non-Gaussian random fields, analytic formula with effects of the weakly non-linear evolution are also derived (Matsubara 
2003). Using the Edgeworth-like expansion, the formulae is given by 
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where 5*"' (a = 0, 1 , 2) are the skewness parameters which are defined by 
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The above formula is a general one for weakly non-Gaussian fields. 

In second-order perturbation theory, the skewness parameters of a smoothed density field, convolved with a Gaussian smoothing 
kernel are given by 
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where the factor S'^fiR) is defined by 
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and is a smoothing length, /iy(z) is the modified Bessel function. In linear theory, the variances CTq and erf are given by 

■ k^dk . 



m+l/2(^l'2), 
2 „„j „2 
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with 7 = 0,1. In this paper we only consider the Gaussian smoothing kernel Wr{x) = tt ^^^R ^ exp(-x^/ 2R^) and thus the window 

function is given by W(/t/;) = e-(**) which is a Fourier transform of the smoothing kernel. The notation Punik) indicates the linear 
power spectrum. The constant E depends on the cosmological parameters and is approximately equal to 3/7 (Matsubara 2003) for 
many sensible models. The values of the skewness parameters are independent on the amplitude of the linear power spectrum in a 
lowest-order approximation. The linear power spectrum is given by the adiabatic CDM model with Harrison-Zel'dovich spectrum 
(Bardeen et al. 1986). Neglecting other components, such as baryons and neutrinos, the skewness parameters S^"\a = 0, 1,2) are 
functions of only a combination TR where F is the shape parameter of the CDM transfer function. The shape parameter F = 0.2 
is adopted in our A^-body simulation of this paper. In Tabled the values of do, cti in linear theory, and the skewness parameters in 
second-order perturbation theory are shown for various smoothing lengths, assuming the same cosmological parameters as that 
of the simulation. 

In Figure^ we plot the functions of equation (|9}. The skewness parameters in this Figure are calculated by assuming TR = 4.0. 
The curves of <to = 0.1,0.2,0.3,0.4 are depicted. The curves with the label ctq = 0.0 correspond to the Gaussian limit (Tq 0. 
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TABLE 1 

The analytical values of o-q, Ro-i/ao and the skewness parameters S*"' calculated in the CDM model. The shape parameter r = 0.2 

AND an approximation £■ = 3/7 ARE ADOPTED. THE POWER SPECTRUM IS NORMALIZED BY =0.9. 
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Fig. 1 . — The MFs for Gaussian and weakly non-Gaussian random fields. The solid lines with label uo = correspond to a Gaussian field. Other curves with 
(To = 0.1,0.2,0.3,0.4 correspond to the weakly non-Gaussian field. 



2.3. Conversion of the Threshold Density by Volume Fractions 
Although we have used the threshold density f = ^th/co so far, we can use another definition of the threshold Vf (GMD) by the 



volume fraction of the high-density regions: 
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(18) 



In random Gaussian fields, the two definitions of v and Vf are identical. In most of the previous work on the genus analysis, 
the rescaled threshold Vi is commonly employed. This re-scaling removes exactly the effect of the evolution of the one-point 
probabiUty distribution of the density field. Since the two-point characteristics of the density field only affect the MF curves 
through the parameters ctq and cti, the remaining differences must be due to higher-order characteristics. 
Analytic formulae of the MFs by the threshold is given by (Matsubara 2003) 
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The highest-order Hermite polynomial in the weakly non-Gaussian terms vanishes in each case. Moreover, the skewness param- 
eters only appear in the form of differences, 5^°' -5^°' (a = 1,2). Therefore, the resultant equation with i/f is simpler than the 
original form using the threshold v. Because the values of skewness parameters are usually close, the non-Gaussian correction is 
small with the rescaled threshold Vf. 



3. COMPUTATIONAL METHODS 
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TABLE 2 

Model parameters of our W-body simulation. 



Model No. of pailicles Physical box size Qq Aq h ag F 
ACDM 512^^ lQ24r'Mpc 0.30 0.70 0.6667 0.90 0.2 




Fig. 2. — The dependence of Euler characteristic on the number of grid points. The smoothing length is fixed to 20/i" Mpc. One-dimensional grid numbers 
are 16, 32, 64, 128, 256 and 512 as indicated in the figure. The curve does not converge and behave awkwardly when we adopt the grid number of 16 or 32. On 
the other hand, the curves of 128, 256, 512 are difficult to distinguish in the figure, thus the result is converged. 



Our A^-body simulation is similar to the one carried out by Jing, & Suto (2002), but the physical box size 1024/i"'Mpc is much 
larger than that. The number of particles is the same, = 5 12^. This simulation was carried out with the particle-particle-particle- 
mesh (P^M) code on the vector-parallel machine VPP5000 at the National Astronomical Observatory of Japan. 

The assumed cosmological model is a flat CDM model with a cosmological constant (ilo = 0.3, Ao = 0.7). The primordial 
density fluctuations are assumed to be random Gaussian, and the power spectrum is given by the Harrison-Zel'dovich spectrum. 
The linear trans fer functio n for the power spectrum of the dark matter is that for adiabatic CDM fluctuations. We use the fitting 
formula by Bar deen et al.l ( 1986) with the shape parameter F = ^oh = 0.2. The amplitude of the power spectrum is determined by 
(78 = 0.9 at the present time. Table|2]shows the parameters in this simulation. 

In studying MFs, we first generate a continuous density field of galaxies. In order to handle the continuous field, we set the 
grids in the simulation box. Then the field smoothing is applied in Fourier space, since the direct convolution in real space is 
much slower to compute. The larger the grid number is, the more precisely the density field is approximated. However, the large 
number of points requires a long time to compute. We need to find an appropriate number of grid points, keeping the calculation 
as accurate aspossible. 

In Figure |2| the Euler characteristic calculated from the A^-body simulation are plotted. The smoothing length is fixed to 
20/!~'Mpc, and the number of grid is varied from 16^ to 512^'. The values of the Euler characteristic converge when the grid 
number is large as expected. Since 256^ and 512^ grids essentially give the same result, we consider results with 256^ grids are 
accurate enough. The 128^' grids already give reasonable approximation. Similar tendencies are seen in other MFs. We conclude 
that the grid number is sufficient when the spacing of the grids is larger than 7?/2.5. Therefore, as long as R > 10/!~'Mpc, 256^ 
grids in our simulation box (1024/i"'Mpc)-' are enough. 

Next, we explain error estimation methods in calculating the MFs from the data of the A^-body simulation. The comparison 
between the analytic formulae in weakly non-linear regions and the MFs calculated from the A^-body simulation can be properly 
done only if we estimate the correct errors in the simulation. There is a constraint in this work that we have only one set of the 
A^-body simulation data. If there is only one set of data, we can obtain only one set of the MFs. To calculate the errors, first we 
divide the cubic box of the A^-body simulation into eight sub-cubes, and then obtain the values of MFs in each cubes. The errors 
are estimated by their averages and variances divided by a scaUng factor, \/7. 
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TABLE 3 

The values of ag, Rcri/cro AND THE SKEWNESS PARAMETERS S*"' FROM NUMERICAL CALCULATIONS OF EQUATIONS ITol-i72l BY THE W-BODY 

SIMULATION. 
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Fig. 3. — The values of cro, i? X cti/cto and the skewness parameters S*"'. The analytical values are shown with hues, and the values of the W-body simulation 
are shown with symbol. 1 a error bars are also shown for the parameters of the A'-body simulation. 



4. RESULTS 

4.1. The values of the skewness parameters 

Before examining the MFs, we first calculate the skewness parameters from the simulation, since these parameters are essential 
in analytic formulae of the MFs. Table|3land Figure|3]show the standard deviation do of the density distribution of the galaxy 5, 
RcTi I Co, and the skewness parameters of equations ( I10t -( ll2t calculated numerically from the A^-body simulation. We use 256^ 
grids commonly for each smoothing length. In Figure |3] the analytic predictions and numerical values from the simulation of 
these parameters are compared. 

As obviously seen in Figure|3j most of the figures agree within the range of numer ical errors. There have been similar studies 
comparing analytic and numerical values of 5*"' and disagreements are reported (e.g.,'Collev et al. 2000; Hikage. Taruva. & Sut^l 
l2003b) . One of the reason of the disagreements is that those samples are not large enough in volume and hence the calculations 
of the skewness parameters are affected by the cosmic variance. Other reason is that the previous work uses biased sample in the 
simulation. Our comparison is directly made by the dark matter distribution in order to separate the complex biasing effects and 



7 




Fig. 4. — The comparison of the MFs for the threshold v at the smoothing length R=i5h Mpc. Thick solid lines: the MPs calculated from A'-body simulation; 
thin solid lines: 1 a errors; short-dashed lines: analytical curves for Gaussian random field; long-dashed lines: analytical curves for weakly non-linear evolution. 



to concentrate on purely nonlinear evolution. 

There are some disagreements in the values of ctq and cri , which are expected if we take into account the fact that non-linear 
corrections of order CTq should be added to the linear predictions of (Tq. Without the second-order correction, the theoretical values 
of CTo have a tendency to the underestimation. This can be understood by the fact that the correlation function integrated over 
whole space have to be zero while the powers on small scales are enhanced by nonlinear evolution. Consequently, the powers on 
weakly nonlinear scales are suppressed and the linear predictions of ctq overestimates the power (Peacock, & Dodds 1994, 1996). 

4.2. The comparison of the Minkowski Functionals between the N-body simulation and the analytic formulae 

In this section, we compare the MFs calculated from the A^-body simulation with the analytic formulae. The calculation method 
of t he M Fs of the A^-body simulation is based on Crofton's formula (Crofton 1868), derived from integral geometry (for details, 
seel Schmalzing & Buchert 1997) . Figures|3|5]and|S]show the four MFs Vkiy) (k = 0, 1,2,3). Each MF is calculated adopting 
smoothing lengths of 15, 25, and 50/i~'Mpc. The analytic formulae of the Gaussian random fields are plotted with short-dashed 
lines, and those of the weakly non-linear fields with long-dashed lines. The amplitudes of the analytic curves are estimated by 
calculating ctq and cti directly from the A^-body simulation. 

Figure0shows the differences between the MFs of the A^-body simulation and that of the analytic formulae. In the plots, the 
differences divided by the 1 a errors of A^-body data are shown. Thus if the curves are in the range of (- 1 , + 1 ), the two curves are 
consistent with Icr significance. As expected, the MFs calculated from the A^-body simulation agree with the analytic formulae 
for large smoothing lengths. 

For quantitative comparisons, we calculate the mean differences between normalized MF values of the analytic prediction and 
that of the simulation: 
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^E[n^'^(-.)-vfv,-) 



(20) 
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where superscript (A) stands for "analytic" and (S) stands for "simulation", and Vi:(max) is the maximum value of a MF Vi 
evaluated by the Gaussian formula, i.e., Vit(max) = V()(i^ = -oo) = 1, Vi(i' = Q), Yi(y= 1) and V3(iy = 0) for fc = 0, 1, 2, 3, respectively. 
The set of various thresholds Vi are given by choosing A^ = 121 equally spaced points with interval = 0.05 in the range 
-2<v< +4. 

When the mean differences defined above are within errors of the numerical simulation, the indicated differences are not 
distinguished from real differences between analytic predictions and true values of MFs. Therefore, we also define the mean 
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Fig. 6. — Same as figure|l] but for the smoothing length R = 50h 'Mpc. 
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Fig. 7. — The ratio of difference between the MFs of W-body simulation and the analytical values to the en'or at the smoothing lengths R =15, 25 and 
50/i^'Mpc. Upper left panel: the volume fraction; upper right panel: the surface area; lower left panel: the integrated mean curvature; lower right panel: the 
Euler characteristic. Xp/d.o.f. and Xg/d.o.f. calculated in the range — 2 < < -1-4 are also shown. 



errors of the normalized MFs in simulation data. 



y,(max)^ A?^ 



(21) 



and the chi-square per degree of freedom: 



1 J!, \vl^\yi)-yf\^i) 



i=\ 



(22) 



where /Sy^\vi) indicate the Ict errors in the simulation data. When Dk < Eu, the indicated differences are within the errors 
originated from a numerical resolution of the simulation, and are not considered as real differences. In fact, this criterion is 
consistent with a usual chi-square criterion < 1. When > 1, the differences are interpreted as deviations of the analytic 
predictions from true values of MFs. Thus accuracies of the analytic predictions are indicated by the quantity when > 1. 
On the other hand, when ^ 1, the statistic Dk represent the upper limits of the accuracies. 

In Table|3 the above statistics are shown for all MFs. In the Table, of^'^^ and D^''"' indicate the equation \2Q\ for 2nd-order and 
linear perturbation theory, respectively. Likewise, Xm xL indicate the equation (I22> for 2nd-order and linear perturbation 
theory, respectively. 



Overall, since Z)*^"''' are smaller than Dj,""' in each case, the second-order predictions show better agreement with simulation 
data. For example, a prediction of the second-order theory of V3 agrees with the simulation data by 4% accuracy level, while that 
of the linear theory agrees by 10% level. In other words, the chi-square statistics of the second-order theory are always smaller 
than that of the linear theory. Since the linear theory is less accurate approximation than the second-order theory, this is a natural 
consequence. When the smoothing lengths are large, the differences between second-order predictions and the simulation data 
are not distinguishable. There are two reasons for this tendency. First, the analytic predictions are more accurate on large scales. 
Second, the simulation data suffers more errors by finite volume effects. 

On smaller scales, R ^ 10/;~'Mpc, the analytic predictions of the perturbation theories are less accurate because of the nonlinear 
evolutions. Even in this regime, the second-order theory is more accurate than the linear theory. Therefore, the predictions of the 
second-order theory is always better than that of the linear theory in the regime we investigate, i.e., 10/;~'Mpc < R < 50/!~'Mpc. 

At this point, one should note that there are correlations between the adjacent thresholds i^, and vi+i . If the value of the MFs 
increases (decreases) at some threshold density f,, there is higher possibiUty that the value also increases (decreases) at the next 



-,(lin) 
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TABLE 4 

Values of Dt of equation t2Ql . Ek of equation i2Tl . and of equation t22l . 









Smoothing Length R(h ' 


Mpc) 








10 


15 


20 


25 


40 


50 




„(2nd) 
k 


0.27% 


0.050% 


0.016% 


0.010% 


0.013% 


0.017% 




r^(!inear) 

k 


1.58% 


0.70% 


0.38% 


0.24% 


0.14% 


0.15% 





Ek 


0.003% 


0.003% 


0.005% 


0.010% 


0.030% 


0.056% 




^2nd 


3631 


34.76 


5.950 


1.371 


0.314 


0.366 






13317 


278.9 


66.29 


19.81 


3.608 


2.850 




k 


7.6% 


3.0% 


1.8% 


1.6% 


2.0% 


1.8% 




T-v(linear) 

k 


12% 


6.9% 


5.2% 


5.2% 


5.1% 


4.4% 


1 


Ek 


0.75% 


0.73% 


0.99% 


1.9% 


5.5% 


6.3% 






723.9 


31.22 


5.284 


1.088 


0.136 


0.163 






6678 


180.2 


60.63 


15.60 


1.610 


1.072 




^k 


14% 


6.7% 


4.0% 


3.1% 


3.8% 


6.1% 






25% 


18% 


13% 


11% 


9.1% 


9.6% 


2 


Ek 


1.6% 


2.1% 


3.0% 


4.3% 


9.2% 


15% 




X2nd 


148.1 


24.71 


3.547 


0.775 


0.204 


0.225 






3146 


128.2 


49.34 


12.10 


1.848 


1.066 






14% 


7.1% 


4.5% 


3.9% 


5.4% 


7.7% 




^k 


18% 


13% 


10% 


8.9% 


9.3% 


11% 


3 


Ek 


2.3% 


3.4% 


4.8% 


6.1% 


1.0% 


14% 




^nd 


265.2 


14.70 


2.164 


0.601 


0.384 


0.297 






841.0 


41.67 


13.06 


5.043 


1.513 


0.784 



TABLE 5 

Same as Table|3] but for the rescaled threshold uf. 



Smoothing length Mpc) 

k 10 15 20 25 40 50 

1.4% 0.45% 0.35% 0.64% 1.5% 1.3% 

^tfinear) j 2% 0.42% 0.39% 0.69% 2.1% 2.4% 

1 Ek 0.75% 0.73% 0.99% 1.9% 5.5% 6.3% 
xLd 4.441 0.457 0.121 0.115 0.083 0.070 
xl 1.917 0.377 0.140 0.126 0.113 0.091 

£)tf"d) 3 1.4% 1.9% 3.4% 5.5% 

^(fmear) j j 1.8% 2.5% 4.3% 6.1% 

2 Et 1.6% 2.1% 3.0% 4.3% 9.2% 1.5% 
xLd 7.397 0.687 0.248 0.206 0.221 0.196 
xL 0.961 0.311 0.303 0.298 0.293 0.285 

of""' 3.2% 2.1% 2.3% 2.6% 5.3% 7.5% 

^Oinear) j^^^ 2.2% 3.1% 3.5% 6.5% 88% 

3 Ek 2.4% 3.6% 5.0% 6.2% 11% 15% 
xLd 6.115 0.758 0.280 0.228 0.519 0.340 
Xv 0.351 0.464 0.408 0.330 0.599 0.431 



threshold density because of the correlation. This is particularly prominent in the Euler characteristic V3. Since the chi-square 
analysis cannot take such correlations into account, the discrepancies between the theory and the data might be over emphasized. 

4.3. The comparison of the MFs in terms of the rescaled thresholds 

In the above comparisons, we use the threshold v = i5/(To- In the following, we use the rescaled threshold Vf and make the 
same comparisons. Figuresl8ll9land[TolcorresDond to Figures|4]|5]and|6l respectively. The differences of the MFs between the 
analytic formulae and the simulation are also shown in FigureFTTI 

From these Figures, we find that the differences in terms of V{ are smaller than those in terms of v for all MFs. Mathematically, 
one of the reason for this is that the skewness parameters only appear as combinations of the form 5*"' - 5*"^ in the analytic 
formulae (eq. 1191 ), and the three values of skewness parameters in CDM cosmological model are close. Since 5*'* -5*^ is 
almost zero in broad range of the smoothing length, the deviations from analytic formulae are mainly depend on the value of 

5(2)_5((« 

The differences D^, errors E^, and chi-square per d.o.f. are also calculated in terms of Vf. These statistics are shown in 
Table |5l It is seen that the differences Dk are within several percent in almost all MFs. Most of the are less than 1 and the 
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Fig. 8. — Same as Figure|3 but for the threshold density ff rescaled by volume fraction. 
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Fig. 10. — Same as Figure |4| but for the smoothing length R = 50h 'Mpc the threshold density Uf rescaled by volume fraction. 
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Fig. 1 1 . — Same as Figure 171 but for the threshold density Vf rescaled by volume fraction. 
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prediction of the analytic formulae are not distinguishable with the simulation results. Interestingly, the second-order predictions 
of Vi, V2 and V3 for R = 10/!~'Mpc are less accurate than the linear predictions. On this scale, both the linear theory and the 
second-order theory can be inaccurate. A possible interpretation is that the strongly nonlinear dynamics affect the MFs to cancel 
the second-order correction term in the analytic formulae. 

As a result, the linear theory and the second-order theory are difficult to be distinguished when rescaled threshold i^f is used. 
This is a good news for testing the primordial non-Gaussianity in the density fluctuations since the nonlinear dynamics in course 
of the density evolution does not significantly alter the shape of MFs in Gaussian fields. 

5. CONCLUSIONS AND DISCUSSION 

We analyze the Minkowski functionals using a large A^-body simulation of a ACDM model with a box size of (1024r'Mpc)^ 
and canonical cosmological parameters, ilo = 0.3, Ao = 0.7, h = 0.6667, F = 0.2. The detailed analysis of the MFs with such a 
large A^-body simulation is unprecedented. The validity levels and regions of the analytic formulae of the MFs are studied in 
detail. We focus on the transition scales from linear to nonlinear evolution, 10-50/i~'Mpc, and calculate the typical differences 
of the MFs between the analytic formulae and the A^-body simulation. 

The variance parameters ctq, cti , and the skewness parameters S^^^a = 0, 1 , 2) are calculated from the simulation and compared 
with analytic predictions. It is found that S^"\a = 0, 1,2) agree with the analytic values within the simulation errors in a large 
range of smoothing lengths. In the previous work of such comparisons, physical sizes of the simulations are not large enough 
as in this work. Thus the skewness parameters from the simulation suffer from the cosmic variance and disagree with analytic 
predictions. It is a new result that the skewness parameters calculated from such a large A^-body simulation as (1024/!"'Mpc)^ 
agree with the analytic predictions. On the other hand, the variance ctq calculated from the A^-body simulation is slightly smaller 
than the linear prediction naturally expected in weakly non-linear regime. 

In this paper we use two definitions of the threshold density, i.e., ly defined by the value of the density field, and i/f defined by 
the volume fraction. Nonlinear effects on MFs are stronger against the threshold ly than i/f . In most of the previous studies on the 
topology of the large-scale structure, the threshold density i^f has been adopted. The reason for this choice is that shape of the 
genus curve is empirically known not to be much affected by the nonlinear evolution. This was also analytically shown in weakly 
nonlinear regime. In this paper, it is shown that the deviations of the MFs on nonlinear scales are indistinguishable with Gaussian 
predictions even though we use a large A^-body sample. To detect the effects of the weakly nonlinear evolution in i^f case, a larger 
sample with less cosmic variance will be needed. Sampling the galaxies in over l(/i~'Gpc)'' comoving volume is within reach 
of the future surveys. For example, in the SDSS Luminous Red Galaxy survey (Eisenstein et al. 2001), the large-scale structure 
up to z < 0.5 is being probed. The covered volume is approximately l(/i"'Gpc)'' in this survey. Therefore, our estimates of the 
errors of MFs in the simulation roughly correspond to the errors in this kind of surveys. If we have still larger surveys, e.g., a 
z < 1 survey with a large sky coverage, the MFs can be measured in exquisite detail. When some non-Gaussian behavior in MFs 
is ovservationally found, it is important to distinguish the non-Gaussianity of the primordial density field from that caused by 
non-linear evolution. This work gives useful information in this respect. 

In this work, we use the density distribution in real space. In redshift surveys, the observed distributions of galaxies are in 
redshift space. The shape of MFs as functions of the threshold is unaffected by peculiar velocity effect in redshift space in linear 
theory (Matsubara 1996). Analytic formula of MFs in redshift space with nonlinear effect has not been derived yet. In future 
work, it will be interesting to investigate weakly nonlinear effect on MFs in redshift space, both analytically and numerically. 

TM acknowledges support from the MEXT, Grant-in- Aid for Encouragement of Young Scientists, 15740151, 2003. YPJ is 
supported in part by NKBRSF(G 19990754) and NSFC. 
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